Observation of pseudogap behavior in a strongly interacting Fermi gas 
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Ultracold atomic Fermi gases present an oppor- 
tunity to study strongly interacting Fermi sys- 
tems in a controlled and uncomplicated setting. 
The ability to tune attractive interactions has led 
to the discovery of superfluidity in these systems 
with an extremely high transition temperature 
p], I2I], near ^ = 0.2. This superfluidity is the 
electrically neutral analog of superconductivity; 
however, superfluidity in atomic Fermi gases oc- 
curs in the limit of strong interactions and defies 
a conventional BCS description. For these strong 
interactions, it is predicted that the onset of pair- 
ing and superfluidity can occur at different tem- 
peratures [3-5]. This gives rise to a pseudogap 
region where, for a range of temperatures, the 
system retains some of the characteristics of the 
superfluid phase, such as a BCS-like dispersion 
and a partially gapped density of states, but does 
not exhibit superfluidity. By making two inde- 
pendent measurements: the direct observation of 
pair condensation in momentum space and a mea- 
surement of the single-particle spectral function 
using an analog to photoemission spectroscopy 
jsl, we directly probe the pseudogap phase. Our 
measurements reveal a BCS-like dispersion with 
back-bending near the Fermi wave vector kp that 
persists well above the transition temperature for 
pair condensation. 

In conventional superconductors, fermion pairs and su- 
perconductivity appear simultaneously at Tc. The single- 
particle, or fermionic, excitation spectrum of a conven- 
tional superconductor follows a BCS dispersion given by 



E,=fi± v/(efe - + A2, 



(1) 



where Ck = Ti^k!^ /2m^ ji is the chemical potential, and A 
is the superfluid order parameter. The lower branch of 
the dispersion (minus sign in Eqn. 1), which is the occu- 
pied one at low temperature, has a positive slope at low 
momentum and then turns around and has a negative 
slope at high momentum. This "back-bending" behavior 
arises because of the excitation gap and is a characteristic 
signature of superconductivity. In unconventional super- 
conductors, such as high Tc superconductors, this back- 
bending in the dispersion has been observed both below 
and, remarkably, above Tc [H. The observation of 
back-bending above Tc represents a dramatic departure 
from conventional BCS theory, which predicts that the 



normal state above Tc is a Fermi liquid with a monoton- 
ically increasing single-branch dispersion and a density 
of states that is smooth through the Fermi surface. The 
departure from a conventional BCS description above Tc 
in the form of a gapped excitation spectra is the essence 
of the pseudogap. 

A satisfactory explanation of the origin and nature of 
the observed pseudogap phase in high Tc superconduc- 
tors has remained elusive due to the complexity of the 
materials. In contrast, ultracold atomic gases are rel- 
atively free of complexity, for example, having no un- 
derlying lattice structure, impurities, or domain bound- 
aries. Moreover, the interactions responsible for pairing 
and superfluidity in ultracold atom gases are well un- 
derstood at the few-body level. As a result, these sys- 
tems are ideally suited for investigating the prediction 
of a pseudogap phase due to pre-formed pairs. There 
is much scientific literature on the topic of the pseudo- 
gap in strongly interacting atom gases with wide-ranging 
viewpoints and conclusions [4, .9.-1 8j , including a recent 
article that predicts no pseudogap phase at all [19]. In 
some theories, the pseudogap phase is predicted to have 
a BCS-like dispersion (Eqn. 1), but where A is no longer 
the superfluid order parameter but instead corresponds 
to an excitation gap due to the formation of incoherent 
pairs 0, lil-llil. However, for the atomic gases, there is 
not yet experimental data to establish the existence of 
a pseudogap phase and confirm its properties. Radio 
frequency (rf) spectroscopy experiments that probe exci- 
tations have been performed above and below the critical 
temperature [2o|, [2l| , but their interpretation relies on as- 
suming a specific dispersion relation and therefore they 
cannot be used to distinguish between a pseudogap and 
a normal phase [13, [ill . 

The question is then, do we have the necessary mea- 
surement tools to look for pseudogap physics in the neu- 
tral atom gas system? To probe the defining properties 
of a pseudogap regime one needs both a measurement 
of the transition temperature as well as a probe of the 
single-particle excitation spectra. In the atomic gas sys- 
tem, the onset of the superfluid phase is clearly detected 
through the observation of momentum-space condensa- 
tion of atom pairs [22]. To probe the single-particle ex- 
citation spectra, we use a technique recently developed 
for atoms that uses momentum-resolved rf spectroscopy 
to realize an analog of photoemission spectroscopy [6|. 
Using these two measurements: the direct observation 
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FIG. 1: Photoemission spectra throughout the pseudogap regime. Spectra are shown for Fermi gases at four different 
temperatures, each with an interaction strength characterized by (kFo)~^ ~ 0.15. The intensity plots show the fraction of 
out-coupled atoms as a function of their single-particle energy (normalized to Ef) and momentum (normalized to /cf), where 
E — corresponds to a non-interacting particle at rest. The spectra are normalized so that integrating them over momentum 
and energy gives unity. White dots indicate the centers extracted from gaussian fits to individual energy distribution curves 
(traces through the data at fixed momentum). The black curve is the quadratic dispersion expected for a free particle, a At 
T — 0.74 Tc, we observe a BCS-like dispersion with back-bending, consistent with previous measurements £]. The white curve 
is a fit to a BCS-like dispersion, Eqn. 1. b,c At T = 1.24 Tc and T = 1.47 Tc, respectively, the dispersion with back-bending 
persists even though there is no longer any superfluidity, d At T = 2.06 Tc, the dispersion does not display back-bending in 
the range of < /c < 1.5 /cf- In all the plots there is a negative dispersion for k/kp > 1.5. We attribute this weak feature (note 
the log scale) to a l/k"^ tail and not to the gap. 



of pair condensation to determine Tc, and, momentum- 
resolved rf spectroscopy to probe the pairing gap, we can 
now explore the issue of the pseudogap in atomic systems. 
In this paper, we report that a BCS-like dispersion, with 
back-bending near kp^ indeed persists even for tempera- 
tures substantially above the measured critical tempera- 
ture for superfluidity. For the atomic gas system, which 
is clean and simple in comparison to high Tc materials, 
this result demonstrates the existence of a pseudogap re- 
gion where incoherent pairs of correlated fermions exist 
above Tc. 

To perform these experiments, we cool a gas of 
fermionic ^^K atoms to quantum degeneracy in a far de- 
tuned optical dipole trap as described in previous work 
[1]. We obtain a 50/50 mixture of atoms in two spin 
states, namely the |/,m/) = |9/2, -9/2) and |9/2, -7/2) 
states, where / is the total atomic spin and m/ is the pro- 
jection along the magnetic-field axis. Our final stage of 
evaporation occurs at a magnetic field of 203.5 G, where 
the 5- wave scattering length that characterizes the inter- 
actions between atoms in the 1 9/2, —9/2) and 1 9/2, —7/2) 
states is approximately 800 ao, where ao is the Bohr ra- 
dius. At the end of the evaporation we increase the in- 
teractions adiabatically with a slow magnetic-field ramp 
to a Feshbach scattering resonance. 

To vary the temperature of the atom cloud, we ei- 
ther truncate the evaporation or parametrically heat the 
cloud by modulating the optical dipole trap strength at 
twice the trapping frequency. To determine the temper- 



ature of the Fermi gas we expand the weakly interact- 
ing gas and fit the momentum distribution to the ex- 
pected 2D distribution [1] and extract (T/Tf)o, where 
the subscript (o) indicates a measurement made in the 
weakly interacting regime, before ramping the magnetic 
field to the Feshbach resonance. For the data presented 
here, we obtain clouds at final temperatures ranging from 
{T/Tf)o = 0.12 to 0.43 with AT = 1 x 10^ to 1.8 x 10^ 
atoms per spin state. The trap frequencies vary depend- 
ing on the final intensity of the optical trap and range 
from 180 to 320 Hz in the radial direction and 18 to 27 Hz 
in the axial direction. Correspondingly, the Fermi energy, 
Ef^ ranges from h-S kHz to /i-13 kHz, where h is Planck's 
constant. The Fermi energy is obtained from N and the 
geometric mean trap frequency, as Ep = hiy{6N)^^^. 
We define the Fermi wave vector as /cf = ^^^^^ and 
the Fermi temperature diS Tp = Ep/ks^ where ks is the 
Boltzmann constant. It is important to note that the 
trapped gas has a spatially inhomogeneous density, and 
one can define a local Fermi energy, and corresponding 
local Fermi wave vector, that vary across the cloud. 

Momentum-resolved rf spectroscopy realizes photoe- 
mission spectroscopy for strongly interacting atoms [6[, 
much like angle-resolved photoemission spectroscopy 
(ARPES) for strongly correlated electron systems. In 
this spectroscopy, an rf photon flips the spin of an atom 
to a third hyperfine spin state and then the spin-flipped 
atoms are counted as function of their momentum. As in 
ARPES [7], conservation of energy and momentum are 
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FIG. 2: Energy Distribution Curves (EDCs) EDCs are 
obtained by taking vertical traces at fixed k through the pho- 
toemission spectra shown in Fig. 1. We show EDCs between 
= 0.1 (top) and = 1.4 (bottom) for the four data sets 
with T/Tc labeled above each figure. Each plotted EDC is 
an average of EDCs over a range of approximately O.IS/cf- 



The EDC at 



1.0 is shown in bold. Black dots indi- 



cate the centers of the gaussian fits to the EDCs. Each EDC 
is normalized to have an area of unity. Vertical dotted lines 
are placed at the local Ef that corresponds to the estimated 
average density of the gas. 



used to extract the energy and momentum of the fermion 
(which in our case is an entire atom) in the strongly cor- 
related system. A key feature of this measurement is 
that the spin-flipped atoms are "ejected" from the system 
in the sense that they have only very weak interactions 
with the other atoms. This means that the spin-flipped 
atoms have the usual free-particle dispersion, and more- 
over, their momentum distribution can be measured us- 
ing time-of-flight absorption imaging with no significant 
effects of interactions or collisions on the ballistic expan- 
sion. This technique was recently applied to a gas just 
below Tc and revealed a BCS-like back-bending disper- 
sion characteristic of an excitation gap [6]. 

To perform the photoemission experiments on atoms, 
we turn on a short rf pulse to transfer atoms from 
the |9/2,— 7/2) state to the unoccupied and weakly- 
interacting 1 9/2, —5/2) state. We then immediately turn 
off the trap and state-selectively image the out-coupled 
atoms on a CCD camera after time-of-flight expansion. 
The rf pulse is kept much shorter than a trap period 
to ensure that the momentum of the out-coupled atoms 
does not change. The length of the rf pulse limits our 
energy resolution to approximately 0.2Ef' As described 
in our previous work, the intensity of atoms out-coupled 
as a function of momentum for each rf frequency can be 
used to reconstruct the occupied single-particle states ^ . 
With this information, one can determine the occupied 
part of the Fermi spectral function and probe the energy 
dispersion. It is important to note that unlike ARPES 



experiments in condensed matter physics the value of the 
chemical potential is not determined in this experiment. 
Rather, in our plots zero energy corresponds to the en- 
ergy of a non- inter acting atom at rest. 

We present our photoemission spectroscopy data 
studying the pseudogap of a strongly interacting Fermi 
gas in figures 1 and 2. The dimensionless parameter 
that characterizes the interaction strength for this data 
is l/kpCL = 0.15(3), where a is the 5- wave scattering 
length. In Fig. 1, we plot the fraction of out-coupled 
atoms as a function of their single-particle energy and 
momentum for temperatures encompassing the pseudo- 
gap regime. In the intensity plots, white dots indicate the 
centers derived from unweighted gaussian fits to each of 
the energy distribution curves, or EDCs, (vertical trace 
at a given wave vector). The energy dispersion mapped 
out with these fits (white dots) can be contrasted to the 
expected free particle dispersion for an ideal Fermi gas 
(black curve). In Fig. 2 we show the same data plot- 
ted as EDCs for wavevectors ranging from k/kp = 0.1 
to k/kp = 1.4. In order to show the evolution of the 
spectral function from below Tc through the pseudo- 
gap regime the data are shown for four temperatures, 
{T/Tf)o = 0.13,0.21,0.25 and 0.35. 

For the data below Tc (Fig. la), we see a smooth 
back-bending that occurs near k = kp- The white curve 
in Fig. la shows a BCS-like dispersion curve, Eqn. 1, 
discussed above; here, we fit to the white dots for mo- 
menta in the range < k < lAkp- While we cannot 
use this fit to extract the gap and chemical potential in 
a model-independent way due to the harmonic trapping 
confinement, the BCS-like fit is consistent with a large 
pairing gap, on order of Ef^ as expected for a Fermi gas 
near the center of the BCS-BEC crossover [H, 13]. 

A striking feature of the data is that the measured 
spectral function evolves smoothly as the temperature is 
increased above Tc. In fact, in all four of our data sets in 
Fig. 1, we observe a weak signal with a strong negative 
dispersion at high momenta. It has been recently pointed 
out that one expects universal behavior at /c >> /ci? for a 
Fermi gas with short-range, or contact, interactions (23| . 
and, moreover, that this will give rise to a weak, neg- 
atively dispersing feature in the Fermi spectral function 
[24] . Recently, we have directly verified this universal be- 
havior with measurements of the momentum distribution 
and found empirically that the expected l//c^ tail occurs 
for k > l.bkp [iH- Therefore, we attribute the negative 
dispersion seen at large k to this universal behavior for 
contact interactions. While the strength of this feature 
should reflect the state of the system, the negative disper- 
sion for k > 1.5/ci? does not, by itself, provide evidence 
of a BCS-pairing gap [22]. 

In the case of a pairing gap, we expect the spectral 
function to exhibit back-bending for k near kp- To avoid 
effects of the universal behavior at large /c, we consider 
the spectral function for k < 1.5/cf- For the three low- 
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est temperatures, we observe a BCS-like dispersion with 
back-bending behavior that occurs near k/kp = 1. We 
interpret this as evidence for the existence of a pseudogap 
regime above Tc comprised of uncondensed pairs in the 
strongly interacting Fermi gas. At our highest temper- 
ature, (T/Tf)o = 0.35, we observe a dispersion that in- 
creases quadraticahy through the region oik/kp = 1 and 
then appears to jump discontinuously near k/kp = 1.5. 
At the same point, the amphtudes of the Gaussian fits 
to the EDCs also drop sharply (see inset to Fig. 3). 
Thus, we conclude that, as the occupation of the posi- 
tively dispersing feature vanishes, the fits jump to a dis- 
tinct, lower energy feature in the spectral function. As 
discussed above, this lower energy feature is consistent 
with the predicted effect of universal behavior at large k 
on the spectral function (23 |. 

In general, for data taken at finite temperature but still 
in the pseudogap regime one might expect to see popula- 
tion in the excited branch of the BCS Bogoliubov disper- 
sion (plus sign in Eqn. 1). Signal in this branch repre- 
sents thermally populated excitations above the pairing 
gap. The data in the region of 1 < < 1.5 are sug- 
gestive of some occupation in this branch. However, the 
limited signal-to-noise ratio makes it difficult to identify 
the excited branch in our data. In addition, the inho- 
mogeneous density of the trapped gas could make the 
observation of two distinct branches more difficult. To 
be conservative, we fit each of the EDCs to a single gaus- 
sian, and we find this to be sufficient to identify back- 
bending. It will be a subject of further research to see 
if the excited branch can be more clearly observed in a 
momentum-resolved atom photoemission measurement. 

In Fig. 3, we directly contrast the dispersions obtained 
for temperatures T /Tc = 0.74, 1.47, and 2.06, shown in 
black, red, and green, respectively. We compare the ex- 
perimental dispersions (circles) to a BCS-BEC crossover 
theory described in Ref. [4]. To compare to the experi- 
mental data, we fit theoretical EDCs to single gaussians 
to extract the centers; the results are shown as lines in 
Fig. 3. The theory incorporates the trapping confine- 
ment as well as the energy resolution due to the finite 
rf pulse duration. The theory, which gives the expected 
k~^ behavior at high k for the momentum distribution, 
agrees qualitatively with the experimental data. Both 
experimental and theoretical dispersions show smooth 
BCS-like dispersions with back-bending near k = kp for 
temperatures up to T/T^ = 1.47. For T/T^ = 2.06, both 
the experimental and theory show a quadratic disper- 
sion before the signal decays around k = l.bkp^ leav- 
ing a much weaker negatively dispersing feature as pre- 
dicted for a normal gas with contact interactions. The 
apparent disagreement between theory and experiment 
at T/Tc = 0.74 can be attributed to a sharp variation of 
the order parameter with temperature close to Tc. In the 
inset of Fig. 3, we show the amplitudes of the gaussian 
fits to the measured EDCs for each of the three temper- 
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FIG. 3: Single-particle dispersion curves. The fits to the 
EDC centers are shown for the three temperatures in Fig. 1 
and Fig. 2 a,c and d, represented by black, red, and green, 
respectively. We observe a BCS-like dispersion, smooth back- 
bending near ^ = 1, for temperatures below and moder- 
ately above Tc. For the highest temperature, we observe a 
quadratic dispersion near ^ = 1 and a sharp discontinuity 

near = 1.5. The lines are theory curves that include effects 
of the harmonic trap and contact interaction, as described in 
the text. Inset We show the amplitudes from the gaussian 
fits to the EDCs for the same experimental data. The fit am- 
plitudes evolve smoothly for the lower temperatures but jump 
discontinuously for the highest temperature gas. 



atures. 

With theoretical calculations for a homogeneous Fermi 
gas, we find that the strongly interacting gas with pre- 
formed pairs and a normal Fermi liquid have distinct 
spectral functions. Namely, the paired state shows a 
smooth avoided crossing (such as described by Eqn. 1) 
while the normal Fermi liquid exhibits a sharp crossing 
leading to a cusp or apparent discontinuity in the occu- 
pied part of the spectral function. The smooth behavior 
in the measured dispersion at the three lower tempera- 
tures, and the sharp jump in the dispersion at large k 
for the highest temperature data are consistent with this 
theoretical picture. 

To determine Tc, we probe pair condensation in our 
atomic Fermi gas following the procedure introduced in 
Ref. [ill. This technique directly probes coherence and 
has been used to map out Tc as a function of tempera- 
ture and interaction strength [22| . The fact that the Bose 
condensation of fermion pairs corresponds to a superfluid 
phase transition was demonstrated unambiguously with 
the observation of a vortex lattice in a rotated Fermi gas 
below Tc [26|. In addition, the accuracy of the conden- 
sate fraction measurements has been investigated both 
theoretically [27|, |28l] and experimentally [29]. In Fig. 4, 
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FIG. 4: Condensate fraction as a function of temper- 
ature. Using time-of- flight expansion at B = 202.1 G where 
our photoemission experiments are performed, we map out 
the condensate fraction. Temperature is measured in the 
weakly interacting regime before the adiabatic ramp to strong 
interactions. We find (^)o = 0.17±0.02. Note that the den- 
sity of the trapped cloud decreases with increasing distance 
from the trap center, and therefore, in a local density picture, 
even at ^ = 0.17 only the part of the gas at the very center 
of the trap is below Tc. 



we show the measured pair condensate fraction as a func- 
tion of the initial temperature of the Fermi gas. As an 
empirical definition of we use the temperature where 
the measured condensate fraction is 1 %. The value of 1 
% is chosen because, when testing our fits with simulated 
data, we find that we cannot differentiate between a Bose 
distribution above Tc and one with a 1 % condensate frac- 
tion. We find Tc = (0.17 ± 0.02) Tp at l/kpa = 0.15(3), 
and using this we report T/Tc for our photoemission spec- 
troscopy data. 

Contrasting the photoemission spectroscopy data with 
this direct measure of the temperature Tc below which 
the system has coherent pairs, we find that BCS-like 
back-bending persists well above Tc in what we iden- 
tify to be the pseudogap phase. Above the superfluid 
transition temperature, these strongly interacting Fermi 
gases are clearly not described by a Fermi liquid disper- 
sion and the existence of many-body pairing well above 
Tc marks a significant departure from conventional BCS 
theory. It is intriguing to note that our measurements 
are qualitatively similar to ARPES results in high Tc su- 
perconductors [8], even though the atomic Fermi gas is 
a much simpler system that does not even have an un- 
derlying lattice structure. However, high Tc materials 
and ultracold atom systems differ substantially and for 
example it may be important to consider that the atomic 
Fermi gas superfluids have a higher ^ compared to high 



Tc superconductors and are more clearly in the region of 
the BCS-BEC crossover. 



METHODS 

Feshbach resonance location 

To create a strongly interacting gas we ramp the mag- 
netic field after evaporation to a value of 202.1 G at a 
inverse ramp rate of 14 ms/G. The data presented here 
is at the same magnetic field as previous "on resonance" 
results presented in Fig. 3b of Ref. [6], where the value of 
a was based on a measurement of the resonance position 
in Ref. [22]. However, from a new measurement based 
on molecule binding energies determined from rf spectra, 
we find the resonance position to be Bq = 202.20(2) G 
and width to he w = 7.1(2) G. With the new resonance 
parameters and B = 202.1 G, we find that the char- 
acteristic dimensionless interaction parameter l/kpa is 
0.15(3). This corresponds to the region of the BCS-BEC 
crossover where the gas is extremely strongly interacting 
and the superfluid gap is expected to be on order of Ep- 
Note that in the absence of many-body physics, the two- 
body prediction of the molecule binding energy at 202.1 
G is 480 Hz, which is less than O.OS^^f- 

Photoemission spectroscopy 
For the work presented here, we have improved the 
signal-to-noise ratio of the photoemission spectra by a 
factor of four compared to our previous measurement 0] . 
Previously, a limitation to the signal-to-noise was related 
to imaging the out-coupled 1 9/2, —5/2) atoms, which lack 
a closed cycling transition for absorption imaging. Now, 
we transfer the out-coupled atoms to the |9/2,— 9/2) 
state with two rf 7r-pulses. Because the number of atoms 
in the |9/2,— 5/2) state is relatively small, this requires 
that we first optically pump the atoms remaining in the 
|9/2,-7/2) and |9/2, -9/2) states to another hyperfine 
manifold. In this way, we can image the out-coupled 
atoms with the cycling transition for the |9/2, —9/2) state 
without contamination from the much larger population 
of atoms that were unaffected by the rf spectroscopy. Be- 
fore constructing the photoemission spectra, we clean up 
the raw images by setting to zero data at large radii where 
the signal drops below technical noise. 

Density inhomogeneity of the trapped gas 
One can define a local Fermi energy, and correspond- 
ing local Fermi wave vector, that vary across the cloud. 
We can estimate average density of the strongly inter- 
acting gas by taking the average density of an ideal 
trapped Fermi gas at a particular ^ and multiplying 

by (^g^)~^/^. Here, is the measured ratio (at finite 

pot pot 

T) of the potential energy of the strongly interacting gas 
to that of a non-interacting gas [30|. For the the data 
shown in Fig. 1 a-d, the local Fermi energy, in units of 
the previously defined Ep^ that corresponds to this aver- 
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age density is 0.81, 0.69, 0.62, and 0.53, respectively. The 
corresponding local Fermi wave vector, in units of kp, is 
0.90,0.83,0.79, and 0.73, respectively. To give a sense of 
the spread in the local Fermi energies, we note that for 
the ideal trapped Fermi gas, the ratio of the local Fermi 
energy at the average density to that at the cloud center 
is approximately 0.6. 
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